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Abstract. We consider semi-discrete first-order finite difi^erence schemes for 
a nonlinear degenerate convection-diffusion equations in one space dimension, 
and prove an error estimate. Precisely, we show that the Lj^^ difference 
between the approximate solution and the unique entropy solution converges 
at a rate 0{Ax^^^^), where Ax is the spatial mesh size. If the diffusion is 
linear, we get the convergence rate 0(Ax^^^), the point being that the O is 
independent of the size of the diffusion. 



In this paper, we consider semi-discrete finite difference schemes for the following 
Cauchy problem 



where Ut = M x (0, T) with T > fixed, u : Ut — > M is the unknown function, / the 
flux function, and A the nonlinear diffusion. Regarding this, the basic assumption 
is that A' > 0, and thus (1.1) is a strongly degenerate parabolic problem. The 
scalar conservation law ut + f{u)x = is a special example of this type of problems. 
Other examples occur in several applications, for instance in porous media flow ['.)] 
and in sedimentation processes [ ^]. 

Since A' may be zero, solutions are not necessary smooth and one must consider 
weak solutions. These are not necessarily uniquely determined by the initial data, 
and in order to get uniqueness, one considers so-called entropy solutions. The 
framework of entropy solutions makes the initial value problem (1.1) well posed, 
for a precise statement, see Section 2. For scalar conservation laws, the entropy 
framework (usually called entropy conditions) was introduced by Kruzkov [24] and 
Vol'pert [30], while for degenerate parabolic equations entropy solution were first 
considered by Vol'pert and Hudajev [31]. Uniqueness of entropy solutions to (1.1) 
was first proved by Carrillo [4], see also Karlsen and Risebro [ " ]. 

For hyperbolic equations, the convergence analysis of difference schemes has a 
long tradition, we mention only a few references. Finite difference schemes have 
been studied by Oleinik [28], Harten et al. [17], Kuznetsov [2(i], Crandall and Majda 
[.^], Osher and Tadmor [29], Cockburn and Gripenberg. [G], Kroner and Rokyta [_ »], 
Eymard et al. [I I], Noelle [32] as well as many others. 
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1. Introduction 



(1.1) 




{x,t) £ Ut, 
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In the last decade, there has been a growing interest in numerical approxima- 
tion of entropy solutions to degenerate parabolic equations. Finite difference and 
finite volume schemes for degenerate equations were analysed by Evje and Karlsen 
[10, 11, 12, l:')] (using upwind difference schemes), Holden et al. [LS, 19] (using op- 
erator splitting methods), Kurganov and Tadmor [- ] (central difference schemes), 
Bouchut et al. [2] (kinetic BGK schemes), Afif and Amaziane [ ] and Ohlberger, 
Gallouet et al. [27, 16, I 'l] (finite volume methods), Cockburn and Shu [ ] (discon- 
tinuous Galerkin methods) and Karlsen and Risebro [22, 2 ] (monotone difference 
schemes) . Many of the above papers show that the approximate solutions converge 
to the unique entropy solution as the discretization parameter vanishes. 

Despite this relatively large body of research, to the best of our knowledge, there 
does not exist a result giving the convergence rate of the approximate solutions to 
degenerate problems. For conservation laws (very degenerate problems), the con- 
vergence rate for monotone methods has long been known to be Ax^^^ [ ], and 
this is also optimal for discontinuous solutions. For non-degenerate problems, the 
solution operator (taking initial data to the corresponding solution) has a strong 
smoothing effect, and truncation analysis applies. Hence difference methods pro- 
duces approximations converging at the formal order of the scheme. However, all 
estimates depend on A' , and are not available if A' is not bounded below by some 
positive number r]. 

Often, the viscous regualarization 

(1.2) u'i + f{u^%^Aiu^U + r,u^,, 

is used to model the behaviour of first order difference schemes. The rationale 
behind this is that first order schemes for (1.1) are formally second order accurate 
for an equation resembling (1.2). If one can prove convergence, or find a convergence 
rate, such that — >■ u, then often analogous arguments will work for appropriate 
difference schemes. If this convergence has a rate, it is expected that the difference 
scheme will have the same rate. In [ ], Evje and Karlsen showed that 

\\u^{;t)~u{;t)\\LHM)=0{v'^'). 

Also, in [ ], Gallouet et al. showed for the boundary value problem corresponding 
to (1.1) that 

\\u^>{;t)~u{;t)\\mM)=0{v'^'). 

To show the same rate for a difference scheme seems remarkably difficult. The 
main result of this paper is that we prove a significantly lower convergence rate for 
a semi-discrete difference approximation u/\x, 

.L-Mt 

/ \u{x,t) -UAxix,t)\ dx <0 (Ax^/^^) , 

where M is a constant larger than |/'|. The O symbol depends on t, L and the 
initial data uq, but not on the discretization parameter Ax. 

The rest of this paper is organized as follows. In Section 2 we make precise 
the definition of a solution of (1.1), and of uak- Then we list a number of useful 
properties of the (unique) weak solution and of the approximation u^x ■ Finally we 
state our main theorem. Section 3 is devoted to its proof, while we test the practical 
convergence properties of the scheme on a numerical example in Section 4. 
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2. Preliminaries 

Independently of the smoothness of the initial data, due to the degeneracy of 
the diffusion, jumps may form in the solution u. Therefore we consider solutions in 
the weak sense, i.e., 

Definition 2.1. Set Ut = {0,T) x R, a function u{t,x) € ((0, T); Li(M)) n 
L°°{Ht) is a weak solution of the initial value problem (1.1) if it satisfies 

D.l A{u) is continuous and A{u)x G L°°(IIt). 

D.2 For all test functions ip G 25(11^) 



(2.1) // uipt + f{u)(px + A{u)ipxxdxdt ^ 0. 

J JHt 

D.3 The initial condition is satisfied in the -sense 

lim / \u{t, x) — uo{x)\ dx = 0. 
Jr 

In view of the existence theory, the condition D.l is natural, and thanks to this 
we can replace (2.1) by 

(2.2) / / uipt + (/(u) - A{u)x) ifx dxdt ^ 0. 

If A is constant on a whole interval, then weak solutions are not uniquely determined 
by their initial data, and one must impose an additional entropy condition to single 
out the physically relevant solution. A weak solution satisfies the entropy condition 
if 

(2.3) Q{u)t + q(u)x + r{u)xx < in V'{Ut), 

for all convex, twice differentiable functions g : M — > M, where q and r are defined 

by 

(f{u) = g'{u)f'{u), and r'{u) = g'{u)A'{u). 

Via a standard limiting argument this implies that (2.3) holds for the Kruzkov 
entropies g{u) — \u — c\ for all constants c. We say that a weak solution satisfying 
the entropy condition is an entropy solution. 
Let the signum function be defined as 

( 1 cr<0, 

sign(a-) = < cr = 0, 
[l cr > 0, 
and its regularized counterpart, sign^, defined as 



Siguier) 



sign(cr) \a\ > e, 
sin ) otherwise, 



where £ > 0. 

We collect some useful information about entropy solutions in the following, for 
a proof see [13]. 

Theorem 2.1. The unique entropy solution u of (1.1) satisfies 

(2.4) // \u-c\pt+&ig^{u~c){f{u)- f{c))ipx + \A{u)- A{c)\ipxxdxdt>Q, 
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for all constants c and all non-negative test functions in I?'(nT)- Furthermore, the 
following limits hold, provided A is strictly increasing. 



(2.5) 



\u -~ c\ <^t+sign(u - c) (/(u) - /(c) - A{u)x) ip^ 



lim / / |A(u)^| sig4 (^(w) - A{c)) ipdtdx, 



(2.6) lim / / {f{u) - f{c)) A{u), sign^ {A{u) - A{c)) ^ dtdx = 0, 

e^O J JtIt 

for all non-negative test functions ip. 

A common method to show existence of an entropy solution is to consider the 
regularized problem 

(2.7) zi? + / (u"), = {A (u^) + rju)^^ , t > 0, u^O, x) = uo{x), 

where ry is some (small) positive number. This equation is not degenerate, and has 
a unique smooth solution for t > 0. The sequence {u''}^^q is compact in L^(IIt), 
and converges to the entropy solution. In [' '-'], it was established that for t < T 

(2.8) h(t,-)-ii"(t,-)llii(M) <CV^, 

where the constant C only depends on /, A and the initial data uq. Of course, 
is also an entropy solution to (2.7), where the diffusion function is given by 

A'^{u) = A{u) + rju. 

Rather than discretizing (1.1), we shall discretize the regularized equation (2.7), 
and let 77 tend to zero in a suitable manner. Due to (2.8), it suffices to compare 

and our approximate solution. To simplify our notation, we therefore, for the 
moment, assume that A is strictly increasing, with A'[u) > 77 > 0. 

We consider a semi-discrete approximation, where space is discrete, but time 
continuous. Let Aa; be some small parameter, and set Xj — jAx and Xj+1/2 = 
(j + l/2)Aa: for j e Z. Set Ij = (Xj_i/2, a;j+i/2]. The discrete derivatives are 
defined by 



■' Ax 



Our scheme is defined by 
(2.9) 



luj{t) + = D-D+Aj, t > 0, 

= :h Ii,^o{x)dx, 



for j E Z. Here is the Engquist-Osher flux and Aj — A (uj). More precisely, 

for a monotone flux /, the generalized upwind scheme of Engquist and Osher is 
defined by 

Fj^i = F{uj,Uj+i) = /+(uj) + f-{uj+i), 

where 

f+(u) = f{0)+ niax(/'(s),0) ds, f- [u) = mm{f {s),0) ds. 
Jo Jo 
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(2.10) 



D-D+Aj, t > 



With this, we can rewrite our scheme (2.9) as 

Uu,{t) + {D-f+{u,)+D+r{u,)) 

for j e Z and Aj — A [uj ) . 

In order to define an approximation on the whole of IIt, we let u/^x be the 
piecewise linear interpolant given by 

UAx{x, t) = Uj{t) + D^Uj{t) {x — Xj) , for x £ [xj,Xj+i], 

and with a slight abuse of notation we define Uj to be the piecewise constant (in x) 
function 

Uj{x,t) = Uj{t) for X e {Xj_i/2,Xj+i/2]- 

We collect some useful results regarding UAa; and the entropy solution u in the next 
lemma. 

Lemma 2.1. Ifu is the unique entropy solution of (1.1) anduj the Junction defined 
by the scheme (2.9). Then the following estimates hold: 



(2.11) 




^)l L=c(jj 


) < 


L°°(R) 




(2.12) 






) < 


"0 LB.y.(R) 




(2.13) 


\\fH;t)) 


— A{u{-, tj)x\\ ^oc(-j{ 


) < 


ll/("o) - A{uo)x 




(2.14) 


\f«,t)) 


— A{u{-,t))x\g Y (^^^ 


) < 


\f{uo) - A{ua)x 


B.V.{R) 


(2.15) 




II"jWIIlo.(k 


) < 


II"j(o)IIl~(b) 




(2.16) 






) < 


K'j(0)LB.y.(R) 




(2.17) 


-^"^+1/2(0 - 




^ < 


\\F^+,/2iO) - D- 


^A{u,{0)) 


(2.18) 


\F] + l/2{t) - 


D+A{u,{t))\^^^^ 






A{u,{0))\ 



iB.V.i 

In addition the initial error is bounded by 

(2.19) ||uA.(-,0)-Mo|lii(K) < CAx, 
for some constant C which depends only on \uQ\g y ^T^y 

For a proof of this, see [! '>]. Note that if uq and A{uo)x (and Uj{0) and 
D^A{uj{Q))) are bounded independently of rj (and Ax), then (2.13) and (2.17) 
imply that 

(2.20) \\D+Aiu,)\\^^^^^ < C and \\Aiu)xh^^^-) < C, 

for some constant C which is independent of Ax and rj. 
Our main result is the following 

Main Theorem. Let u be the unique entropy solution to (1.1) and be as 
defined by (2.10). Choose a constant 

M > max \f'{u)\ , 

I"I<I|moIIi,oo(5) 

and another constant L > MT, where T > 0. Then there exists a constant C , 
independent of Ax, but depending on f, L, T and uq, such that 

/ \u{t,x) -UAx{t,x)\ dx <CAx'^^^'^ fort<T. 

JL+Mt 
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As a by-product of our method of proof we get an improved rate if the diffusion 
is hnear. The significance of this rate is that is independent of the size of the 
diffusion, which in this case is rj. 

Main Corollary. Let u he the unique solution to the viscous regularization 
ut + f{u)x = rjUxx, i > 0, u{x, 0) = uo(a;), 

and let uax &e defined by (2.10) with A{u) — rju. Then there exists a constant C , 
independent of Ax and rj, hut depending on f, L, T and uq, such that 

pL-Mt 

/ \u{t,x) -UAx{t,x)\ dx <CAx'^/^ fort<T. 

JL+Mt 

In the hnear case this is what we expect. In fact, in [■")] Chen and Karlsen 
showed that for a linear flux function / the expected r] independent rate of 1/2 
holds. However, the methods used in [ ] are not easily modified to nonlinear flux 
functions. 

3. Proof of the main theorem 

First of all, for simplicity we will prove the main theorem for /' < 0. Note 
that, in that case = /(uj+i). The general case treatment will be similar 

(see Remark 3.1). The theorem will be proved by a "doubling of the variables" 
argument, but we start not with the entropy condition (2.5), but in the argument 
leading up to this condition. Set 



i^eiu, c) = y sign^ {A{z) - A{c)) dz. 
This is a convex entropy for all constants c. Set u — u(y, s) and rewrite (1.1) as 



Us + {f{u) - f{c)\ = {A{u) - A{c))yy , 



and multiply this with V'e(wj c)(/7 where is a test function with compact support in 
E X (0,T). Remember that A' > 77 > 0, so m is smooth, after a partial integration, 
we arrive at 

V'e(u, c)ips + Qt{u, c)ipy dyds 

sign^ {A{u) - A{c)) A{u)yipy + sign^ {A{u) - A{c)) {A{u)yf ipdyds. 



Where we have used Q'^{u,c) ~ ^p'g{u,c)f'{u). Although tps{u,c) fa \u — c|, is 
not symmetric in u and c. This makes it cumbersome to work with when doubling 
the variables, so we rewrite the above as 

\u - c\ips + Q^{u,c)(py dyds 



sign; {A{u) - A{c)) {A{u)yf ^ + sign, [A{u) - A{c)) A{u)y^y 

+ (|w - c| - ^e{u, c)) dyds. 

In the doubling of variables argument we choose c = u/y^x{x, t), and a test function 
ip{x,y,t,s). Integrating the above over {x,t) € 11^ after an integration by parts, 
we end up with 
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(3.2) / \u-UAx\^s +sign^iA{u) - A{uAx))if{u) ~ f{uAx))^ydX 
Jul 



{A{u) ~ A{uAx)) ({A{ U)y ) -^A{u)yA{uA 

\\^ (ipyy + ipxy) + (V'e "Aa;) - \u ~ UAx\) 

(sigi\(A(z) - A{uAx))){f{z) - f{uAx)) dz I Lpy 



\A{u) - A {uAx)\^ iVyy + ^xy) + (V's "Aa;) - \u ~ UAx\) 

dz 

where dX — dydsdxdt and \a\^ — sign^[z)dz. Here we have used that 

= / (sign^ {A{u) - A{uAx)) A{u)yip) dX 
Jul 



dX, 



n2. 



sign^ {A{u) - A{uAx)) A{u)yA{uAx)xV> - - ^ (wax)!^ "^xy dX, 



and that 



s\g-!i^{A{u) - A{uAx)) A(u)yipydX ^ / {\A{u) - A{uAx)\^)y'PydX 



11% 



\A{u) - A (max)Ij Vvv 



Also 



,{u,UAx) 



sign^(A(z) - A{uAx))-^{f{z) - f{uAx))dz 



" d 
. dz 



(sign^(A(z) - A{uAx))){.f{z) - f{uAx)) dz 



+ sign^(yl(M) - A{uAx)){f{u) - f{uAx))- 

The next goal is to obtain an analogous estimate for the difference approximation 
uax- Set fj{t) — ip{xj,t) and multiply the scheme (2.9) with ^'^{uj,c)Lpj and do a 
summation by parts to get 

^ (wj, c)j + ipjD+Q,{uj,c) 
j 

= - ^ (%-, c) D+A (u,) + Z?+i^; {uj,c) D+A (u,) 

j 

Where we have used the following result: 

0> / ^,{s,c)[f{u^+i)-fis)]ds= i:^(s,c)f\s)ds 



- V'£(Uj)[/("j + l) - f{Uj)]. 

Set (^^^ = (^j for X £ Xj+1/2], integrating the above for t e [0,T] and 

multiplying with Ax, we obtain 



\uax - c\ " + sign, {A (uj) - A{c)) {f{u,) - /(c)) D"^^- dxdt 
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+ / / (V'e iuj,c) - \uAx - c|) dxdt 

sign^ {A {uj) - A{c)) D+ A {uj) D+ip^"" dxdt 



Hi 

> (( sign; {A (0,+i/2) - A{c)) [D+A {u,)] ' ^,+1 dxdt 

J JyIt 

^{signMiz) - A{c))){f{z) - f{c))dz ] D-^^-dxdt. 

This can be rewritten 

\uAx - c\ipt+ sign^ {A{uj) - A{c)) (/ (uj) - f{c)) (p^ dxdt 



> II sign^ {A {uj) - A{c)) D+A {ujjf^ dxdt 

sign^ (A (6*^+1/2) -A{c)) [D+A{uj)f ipdxdt 
+ 11 sign; {A {0J+1/2) ~ A{c)) [D+A (u^-)] ' iVj+i " f) dxdt 

J JUt 

+ 11 sign^ {A {uj) - A{c)) D+A (uj) {D+ip^'' - ip^) dxdt 

+ 11 {\u/^x - C| - ■(/'e {Uj,c)) tpf"^ + |uAa; " c| [kpt - dxdt 
J JyIt 

+ 11 Sign, {A{u,) - A{c)) {f{u,) - f{c)) - D'^f^) dxdt 

( r ^(sign,(A(z) - A(cmf{z) - /(c)) dz] D-^"^- dxdt. 

In order to make this more compatible with the corresponding equahty for the exact 
solution (3.1), we rewrite again to get (recall that Uj denotes the piecewise constant 
function taking the value Uj{t) in the cell (xj_i/2, a;j_|_i/2]), 

(3.3) 

\uAx - c\ipt+ sign, {A{uAx) - A{c)) (/ (uax) - I{c)) ipx dxdt 

(3.4) > // sign^ (A (uAx) - A{c)) A (u Ax) xVx dxdt 

J Jut 

(3.5) +// sign'^ {A (uAx) -A{c)) [A {u Ax) x]\ dxdt 

+ 11 [signe {A (uAx) - A{c)) - sign, {A{uj) - A{c))] 

J JHt 

(3.6) X (/ (uAx) - /(c)) ipx dxdt 

(3.7) + // sign, {A (uj) - A{c))) (/ (uax) - f (uj)) ipx dxdt 

J JUt 

[sign; {A {0J+1/2) ~ A(c)) - sign; {A {uax) ~ A{c))\ 
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(3.8) X [D+ A {uj)]^ ip dxdt 

(3.9) +// sign'^iAiuAcc)-A{c))\{D+A{u,)f-{A{uA.)J^]ipdxdt 



(3.10) + // sig< {A (0,+i/2) - A{c)) [D+A (u,)] ' (^,+i - ^) dxdt 

J JHt 

[sign, {A (uAx) - A{c)) - sign, {A{uj) - A{c))] 

(3.11) X D+A {uj)ipx dxdt 

(3.12) + / / sign, {A {uA^} - A{c)) [D+A (uj) - A (ua^) J dxdt 

(3.13) +11 sign, {A {uj) - A{c))D+A (uj) {D+ip^'' - ip^) dxdt 

(3.14) +// {\uAx-c\-ip,{uj,c))Lpf'' + \uAx~c\{pt-vt'') dxdt 

(3.15) + // sign, {A{u,) - A{c)) {f{u,) - f{c)) - D-^f^) dxdt 

(3.16) + jj^ ^(sign,(A(z) - A(c)))(/(z) - /(c)) dz) i^-^p^^^ dxdt. 

Remark 3.1. In the general case, we can write Fj+^ = f^i^'^j) + f~{'^j+i)- 
Note that {}+)' > 0, {f~Y < 0, and both are Lipschitz continuous. Then if 
we multiply (2.10) by ^'^{uj,c), numerical flux part can be written as D^Qf + 
D+Q~ + "terms having the right sign". Here {Qfy{u,c) — ip'^{u, c){f+y (u) , and 
Qt + Q7 =Qe. 

Let now 

7^l(c) = (3.6) + (3.7) + (3.8) + (3.9), 

7^2(c) = (3.10) + (3.11) + (3.12) + • • • + (3.16). 

Now we choose c = u{y, s) and the same test function as before, and integrate the 
result over (y, s) £ Ht- The result reads 

(3.17) / \uAx-u\ipt + sigi\{A{uAx)- Mu))if{uAx)- f{u))iPxdX 
Jul 



> / si: 



ign^ {A (uAx) - A{u)) ((A {uAx)S - A (uax)^ Mu)v) V dX 

\A{uAx) - A{u)\^ {ipxx +Vxy) dX 



111% 

+ 11 ni{u)+n2{u)dyds. 

Adding this and (3.2), we get 
(3.18) / \uAx - u\{ipt + (fis)) 

+ sign, {A (uAx) - A{u)) (/ (uax) - /(«)) {fx + fy) 
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> 



> 



+ \A{uAx) - A{u)\e ifxx + ^iPxy + 'Pyy) dX 

/ sign; {A (uAx) - A{u)) {A {uax)^ - A{u)yf dX 

+ {^pe{u,UAx)-\u-UAx\)<fsdX 

+ [ (r Y^^^gnAA{z)-A{uAxmm-fiuAx))dz)^ydX 

+ // ni{u)+TZ2{u)dyds 
J Jut 

" d 



dz 



(sign^(A(z) - A{uAx))){f{z) - f{uAx)) dz Lpy dX 



+ {lps{u,UAx)-\u-UAx\)VsdX 

+ 11 ni{u)+n2{u)dyds. 

=■■ I Qi + Q2dX + II ni{u)+'R2{u)dyds 



Now we are going to speciiy a nonnegative test function (f = ip{t, x, s, y) defined 
in X Ht- Let w e Co°(M) be a function satisfying 

supp(w) C [—1, 1], w(it) > 0, / oj{(j) da = 1, 

and define LVr{x) = u){x/r)/r. Furthermore, let h{z) be defined as 

'O, z<-l, 
h{z)=iz + l z£[-l,0], 
1 z>0. 

V 

and set ha{z) = h{az). Let y <The two numbers in (0, T), for any a > define 

J — OO 

^{x,t) = {H^^it -y)- H^^it - r)) (h^{x - Li{t)) - K{x - Lr{t) - ^) 

where the lines Li^j. are given by 

Li{t) = -L + Mt, Lr{t) =L-Mt 

where M and L arc positive numbers, M will be specified below. With < r < 
min {v, T — t} and ag G (0, min {i> — r,T — t — r}) we set 

(3.19) ip{x,t,y,s) = -^{x,t)u}r{x-y)uJro{t- s). 

We note that cj) has compact support and also that we have, 

'Pt + 'Ps = 'bt{x,t)uJr{x - y)u]ro{t - s), 
Vx + 'Py = "^xix, t) UJr{x - y) LOro{t " «), 
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ifxx + 2^Pxy + (Pyy = xx {x , t) UJr {x - y) UJro " s) • 

For the record, we note that 

^t{x,t) - {h'o.{x - U{t)) + h'^{x - LS) - 



"^.{x.t) = x".%)(t) (h'^{x^ U{t)) - h'^[x -Lr{t)--)\, 



(3.20) 



^xx{x,t) = x".V)(t) [h^{x - - h^{x ^ LAt) - -) j . 

We shah let aU the "smaU parameters" a, ao, r, tq, e and Aa; be sufficiently 
small, but fixed. The goal of our manipulations is to obtain an inequality where 
the difference between and u is bounded by some combination of all these 
parameters. 

We shall repeatedly use the fact that 



(3.21) 



for V = u, V ~ u^x, V — f{u), V — A(u)x or v = A{u/^x)x- These estimates follow 
from the basic bounds in Lemma 2.1. Starting the first term on the left of (3.18), 
we write 

\uAx - u\{lPs + (fit) dX < / \uAx{x,t) - u(x,t)\'^tdxdt 



\v{x, t) — v{y, t)\ ujr{x — y) dxdy < Cr 
and // \v{x, s) — v{x,t)\ujro{'t — s) dxds < Cro, 



\u{x, t) — u{x, s)\ \'^tix,t) \ ujrg (t — s) dsdxdt 



\u{x,s) - u{y, s) \ \^t{x,t)\ujroit - s)ujrix - y) dX . 



To estimate /3 and 7 we use 

l^-tl < (wao(i ->^)+ ^coit + M {h'Jx - Li{t)) + h'^{x^ L,,(i) - , 



and that 



]dt 



j uJaoit -v)dt<l and j Mh'^{x - L;,r(i)) 
A typical term in /3 reads 

\u{x, t) — u{x, s) I uJag {t — v)(^ro ~ ^) dxdsdt 



< c. 



< C J J \t — s\ UJro ~~ ^ao ^ ^) dsdt 

<Cro, 
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Hence 

13 < Cro. 

Similarly a typical term in 7 can be estimated 

\u{x, s) — u{y, s)\ ujr{x — y)uJro{t — s)ujao{t ~ dxdydtds < Cr. 
Thus we find that 

(3.22) /3 + 7< C(ro + r). 

To continue the estimate with the first term on the left of (3.18), we split 6 as 
follows 



(5 - 



^ x"j;,)(t)M (^h'^ix - Li{t)) + h'^ix - Lr{t) ^ \uAAx,t)-u{x,t)\ dxdt 



^ \uAx{x,t) - u(x,t)\ {uJaoit - - UJaoit - t)) 



dxdt . 



The term i5i will be balanced against the first order derivative term on the left hand 
side of (3.18). To estimate S2 we set e(x,t) = \uAx{^,t) ^ ^^(2^1^)! and proceed as 
follows 



Hi 

< 



< 



X(L, ,Lr) (^' ^) '^"o ~ ^) ^^^^ 

+ J J x"l,,l,) (2;, |e(a;,t) ~ e(a;,j^) I WQ,o(t- J/) 

X{L,,Lr) ^)e(a;, i^) dx + Cao 
and similarly 

X{Li,L^)ix^t)^i^^*)^ao{t~T)dxdt> J x"L,,L,,)(a;,T)e(x,T) dx-Cao- 
Using this we get the estimate 

h < [ X{Li.L^)i^^'^) \uAx{x,i^) ~ u(x,v)\ dx 



(3.23) 



Now we rewrite the "first derivative term" on the left hand side of (3.18). Doing 
this, we get 

sign^ [A (uAx) - Mu)) if (uAx) - fiu)) (ipx + ^y) dX 

= / sg{x,y,t,s){f{uAxix,t))-'f{u{x,t)))'i'x{x,t)ujrix-y)ujr„i't-s)dX 
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sg(a;, y, t, s) {f{u{x, t)) - f{u{y, s))) t) ujr{x - y) ujro (t - s) dX 



=: ^3 + S4, 



where we have set sg{x, y,t,s) = sign^ {A{uAx{x,t)) — A{u{y,s))). We proceed as 
follows 



l<54|< / ^ |/Kx,t))-/Ky,.))|x(,%)(t) 

ojroit- s)ujr{x-y) (^h'^{x - Li{t)) + h'^{x - Lr{t) - ^) ] dX . 
Each of these two terras are estimated using (3.21) by 

\f{u{x, t)) - f{u{y, s))| x",V) (t) uJro {t - s) uJrix - y) h'^ix - Li{t)) dX 



< 



t)) - fiuix, s))\ x^,V) it) ujr„(t - s) ujrix ~ y)h'^{x - Li{t)) dX 

\J{u{x, s)) - f{u{y, ^roit - s) u:,{x - y)h'^{x - Li{t)) dX 

<C{ro + r), 
and thus 

(3.24) \Si\ < C [ro + r) . 

The terms di + 63 is bounded by choosing M sufficiently large (all functions are 
functions of (x, i)). 

61 + S3 

x":,4i)K{x - Li{t)) i-M \u^x -u\+sg (/(UA.) - fiu))) dxdt 
+ 11 xl:^^^it)h'^ix-Lrit)--){~M\uA-.^u\-sg{f{u^x)~ fiu))) dxdt. 

Choosing M larger than the Lipschitz norm of / implies that 

(3.25) Si+S3< 0. 
Turning to the "second derivative term" (3.18), we get 



/ \A{uAa^)~ A{u)\^^^ 

Jul. 



. UJra (^ ~ s) w,. (x — y) dX 



< I \A{uAx)-A{u)\(h':^{x~Li{t)) + hl{x~Lr{t)~-) 
Jul \ a 

LUro [t ~ s) UJr{x — y) dX 

<Ca {\uAx{x,t)\ + \u{y,s)\) uJroit- s)uJrix~ y)dX 
Jnl 

< Ca, 

for some constant C which is independent of the small parameters. 
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Collecting (3.22), (3.23), (3.24), (3.25) and the above inequality we see that 

■xf.„.,,(..,.)|..,(....)-.fe.)i.. 



< 



(3.26) 



+ C (ro + r + ao + a) 



Qi + Q2dX 



TZi{u) + TZ2{u) dyds 



In order to estimate the integral involving Q2 we first observe that since A' {u) > rj, 
we get 

|^,(a,6)-|a-6|| < P \sign^{A{z) - A{b)) - sign{A{z) - A{b))\ dz 

Jb 



< 



1 



V JA{b) 



sm 



7r/3 
2e 



d/3 



< C- 



Using this we find that 



(3.27) 



dX < C- 



T i-L-{M/a)t 



dxdt 



foV Jo J-L + {M/a)t 



< C- 



L 



1 



where C is independent of the small parameters. To estimate the integral of Qi we 
use the Lipschitz continuity of / (A^^). Note that / (A^^) is Lipschitz continuous 
with Lipschitz constant where M is the Lipschitz constant for /. 



(3.28) 



(sign^ {A{z) ~ A{uAx))) 



\ifiz) - f iuA.))\ dz 



dz 

A{u) 

Sign; (r - A (UA.)) \ (/ {A-\r)) - f (ua.)) \ dr 

min{A{u) ,A{u/^^)+e) 

(/ {A-'(r))-f{uA.)) \ dr 



min(A(ti),A(uAa;) — e) 



< 



Using this we find that 







/ QidX 






V Jn 



dX 
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Ql + Q2 dX 



Therefore 
(3.30) 

Now we claim that 

// Uiiu) + n2iu) dyds 

(3.31) 



Ce / L 1 

< — L+ — + — 
77 \ vq ra 



< c 


Ax 







1 1 1111 

-IT + ^ + ^ + \ h- 

^e-^ er;^ tq sr r 



1 1 



1 



1 



where C depends on (among other things) L and T, but not on the small parameters 
ao: ct, To, r, rj or £. The proof of this claim is a tedious computation of all the 
terms. We start by the one from (3.6). 

1(3-6)1 L=„ dyds < [ Isign, (A(ua.) - A{uj) - sign, {A{uj) - A{uj)\ 

X |/(wAx) - /(u)| l^ajW^Wro + 'i'Uj'r^rol dX 



< c 



CTAx 

< max \U^^\g y. 



UAx - uA l^ifJ + dxdt 



< 



e te[o,T] 
CAx 
er 



r 
1 

a+ - 
r 



for sufficiently small r and a. Now 

1(3-7)1 |,^„, Ciyrfs < M I \uAx - Uj\ I'i'a^UJrUJro + ^t^^rol dX 



< CTAx max ImasIri// 
te[o,T] 



a + - 
r 



< 



CAa 



To estimate the next term we observe that I sign" (z) I < C/e^ 



1(3.8)1 dyds< 



C 



CT 



^i+i/ 



1/2 - UAx\ dxdt 



< — 7T-Ax max \uax 



te[o,T] 



B.V.( 



< 



CAx 

-2 ■ 



The next term involves D^A{uj) and j4(?iA2;)2:, these can be written 

D+A{uj) = A'{aj+i/2)D+Uj and A(ua^)^ = A'(/3j+i/2)i:'+Mj, 

if a; G [xj,Xj^i). Here aj^x/2 ^nd are between and Uj+i. Furthermore, 

since yl' > r;, we have that 

\D+uA < - \D+A{uj)\ . 
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1(3.9)1 1 dyds <— [[ \D+A{uj) - A{uacc)x\ dxdt 



C 
< — 



V£ JJyit 



\aj+i/2 - Pj+i/2\ dxdt 



< — ^ J J — Uj\ A{uj)\ dxdt 

CT 

CAx 

We continue with 

1(3.10)1 1^^^ dyds 

<C \ \D+A(u,)\ujr^{x-y) 

X |5'(a;j+i,t) - y) - ^(x,t)u:r{x - y)\ dX 

X \^x{z,t)u:r{z ~ y) + '^{z,t)u:'^{z — y)\ dzdxdydt 

< C / 2J |^^^(^tj)| / / \'i'xiz,t)\ + -'^{z,t)dzdxdt 

j Jxj-i Jx ^ 

< CTAx (a+-] sup \\D+A{uj{t))\\ . 

V ^/ te[o,T] ' 

CAx 

r 

for sufficiently small a. With similar arguments we show that 

r'Ar 

1(3.13)1 + 1(3.15)1 dyds<—, 
and 

1(3.11)1 + 1(3.12)1 dyds < ^ + 

The term (3.14) consists of two parts. The first of these 

/ li^eiuAxiU) - |max - "ll'^f'^rf-^ < — + - ) , 

by the same arguments used to show (3.27). For the second part, using the fact 
that {uAx)t G -^^(IR)i we can show that 

\uAx - u\ {ipt - {^j)t) dX 



< 



n|. 



\iuAx)t\ {lp{Xj,t) - 1p{x,t))uJr{x - y)uJro{t " s) 
+ \{uAx)t\ i^^rix -y) - UJr{Xj - y))lp{x ,t)Uro{t - s) 



dX 
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< CAa 



CAx 

r 



There remains the last term, (3.16). To estimate this, we follow the similar 
arguments to the one used in (3.29). The result is that 



1(3.16)1 1 _ dyds<C- 



rrj 



L 



Now we have proved the following lemma: 



Lemma 3.1. Assume that u and u^x o,i^e take values in the interval [^K,K] for 
some positive K. Let M > maXi,g[_x.A'] l/'(^)l- Then iJT>T>v>Q and 
L — Mt > 0, we have 



\uax{x,t) - u{x,t)\ dx 



L+Mt 



< 



(3.32) 



C 



uax{x, v) — u{x, 1^)1 dx 
1 



ro + r + a + Ax 



F I — 



1 



1 



1 1 


1 


1 


— + ^ 


H 

ro 


H 

er 


1 

H h 


-)] 




rr]a 







\uAx{x,t) - v{x,t)\ dx<C[a+^ + ^ 



This follows from (3.26) and (3.31), observing that we can send ao to zero. Now 
we let v{x,t) be the unique entropy solution of (1.1), where A'{v) > 0. We set 
o; — r = tq = r]^^^, and assume that a is sufficiently small, then 

pL-Mt 

(3.33) / 

J -L+Mt 

for some constant C which is independent of the small parameters a, s and Ax. 
This follows from (3.32), (2.8) and (2.19). Setting a = Ax^/", and e = proves 
the main theorem. 

Proof of the main corollary. To prove the corollary, we retrace the proof of the 
main theorem, using that A(u) — rju. We begin by setting 



ipe{u, c) = rj J sign^(z — c) dz. 
In this case, the equation corresponding to (3.2) reads 



(3.34) 



\U - UAx\ + sign [u - UAx) {f{u) - f {uAx)) dX 



— lim 



sign^ [u - UAx) {ul - Uy {uAx)x) f 

-\U- UAx\^ {^yy + <Pxy) 



dX. 



On the other hand, to obtain an expression for uax we can proceed as follows. In 
fact, in discrete set up we have the following inequality 

Tpi;{uj,c){ipj)t + Qe{u.j,c)D^Lpj dxdt 



18 



K.H. KARLSEN, U. KOLEY, AND N. H. RISEBRO 



>ri D^UjD^ipjSign^{uj - c) dyds + V c^a; [sign^ (uai; - c)]{uax)x'P dxdt 



dxdt 



L»+[sign^(wj - c)\D+UjLpj+i - 9^ [sign^ (maj; - c)](wax)xV 
^^JJ i'^Ax)xfxSign^{uAx - c) + sign^(uAa; - c){uAx}lfdxdt 
+ 'n (wAx)K£cSignj(uA2; - c)(p - (pjsign^{uj - c)D^D^Uj dxdt, 



where we have used Q'^ = ^e/' and the following equality: 
0> / ^,(s,c)[/K+i)-/(s)]ds= / ^p,is,c)f'{s)ds 

J Uj J Uj 

By taking limit as e — >■ 0, and choosing c = s) we have 
(3.35) 

\uj - u\ {ipj)t + sign(uj - u){f{uj) - f{u))D^ipj dX 



dX. 



> lim 

e->-0 



{iuAx)l - Uy{uAx)x) Sign^ (uAx -u)ip- \uAx - {ipxx + Vxy) 

Now adding (3.34) and (3.35), using Q[a, b) = sign(a — b){f{a) — f{b)), we get 

I^Ax - U\ {(fit + (fis)) 

+ sign {uAx - u) (/ (uAx) - /(-")) iVx + Vy) 
+ Wax - ""le [Vxx + '^fxy + Vyy) dX 



> i\uj -u\~'\uAx~u\){ipj)tdX + \uAx-u\{ipt~ {ipj)t)dX 
Jul Jul 

+ Q{uj,u){lPx - D^fj) dX + {Q{u,UAx) - Q{uj,u))(pxdX. 

Observe that we can use similar arguments to the ones used in the nonlinear diffu- 
sion case for all the terms on the left side of the above inequality (3.36). Remaining 
all the terms can be estimated in the following manner. We begin with 

{\Uj - u\ - \uAx - u\) {(pj)tdX < / \uAx - Uj\\i^tUJrUJro + 'Ipi^r^rol dX 



CaTAx , , 

< max kiAr D 



y.(R) 



< 



CAx 

ro 



Next, we continue with 

Q{uj,u){(px- D~ipj)dX = - I Q{uj,u)xip - D'^Q{uj,u)(pj dX 



I I Y.I (Ax^5^,^^^^^}D+Q{u„u) - D+Q{u,,u)^,) dX 
Jut Jo Jxj ^ ' 
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CAx 



< 

r 



In a similar way, we also can show that 

(Q(u,-UAa;) - Q{Uj,u)) (fix dX < 



r 



Finally, we end up with 

t) — u{x, t)\ dx 



(3.36) 



L+Mt 



A.T Ax 

< / \uAx{x,v) — u(x,v)\ dx ^ C [ rQ + r + a ^ 1 

r To 



Now setting a = r = tq yields the estimate 

/ Ax 

\uAx{x, t) — v{x, t)\ dx < C \ a -\- 



-L+Mt 



a 



Setting a — Ax^/^ proves the corollary. □ 

4. A NUMERICAL TEST 

In order to test the unlikely optimality of the convergence rate of our main 
theorem, we compute the numerical convergence rate of an example. Consider the 
following initial value problem 

ut — A{u)xx for i > and X G (— 7r/2, tt), A{u) — ^ {max {u,0})'^ 
u(x, 0) = sin(x), X e [— 7r/2, tt], 



(4.1) 



supplemented with the boundary conditions 

dxA{u{t, x)) = for < > and x ~ — 7r/2, x = tt. 

In order to simplify matters, we have chosen an example without the convective 
term, nevertheless a discontinuity will form in the solution due to the degeneracy 
of A. This discontinuity will take the form of a boundary on which u — 0, moving 
to the left. 

We have used the Euler method to integrate the system of ordinary differential 
equations (2.9), resulting in the update formula 

Uj{{n + l)At) = Uj{nAt) + AW^ D+ A.j{nAt). 

For this to be linearly stable, the time-step At must obey the restrictive CFL 
condition At < 0.5 max A'(u)Ax^. 

In Figure 4.1 we show the solution in the {x,t) plane and a snapshot of u at 
t — 1, for an approximation using 400 grid points in the interval (~7r/2,7r). 

Finally we computed approximate errors, defined by 

^qq \\uAx - Mrcflli,! 
IWaxWli 

for a reference solution computed by using our scheme with 4000 grid points. These 
errors were computed at t = 1. The result of this is shown in Table 4.1. Looking 
at this table, it seems the numerical rate of convergence is 1. Thus very far from 
the rather pessimistic lower bound proved in this paper! 
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X 



Figure 4.1. An approximate solution to (4.1) using 400 grid 
points. Left: u in the (x,t) plane for t € [0,4]. Right: an ap- 
proximation to it(l,x) using 25 grid points, a reference solution 
computed using 4000 grid points and the initial data. 



N 


25 


50 


100 


200 


400 


800 


error 
rate 


3.62 


1.55 
1.22 


0.82 
0.92 


0.40 
1.02 


0.18 
1.11 


0.07 
1.42 



Table 4.1. Numerical errors and convergence rate for the intial 
value problem (4.1). N is the number of grid points. 
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